Effects of thermal phase fluctuations in a 2D superconductor: an exact result for the 

spectral function. 
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We consider the single particle spectral function for a two-dimensional clean superconductor in 
a regime of strong critical thermal phase fluctuations. In the limit where the maximum of the 
superconducting gap is much smaller than the Fermi energy we obtain an exact expression for the 
spectral function integrated over the momentum component perpendicular to the Fermi surface. 



PACS numbers: PACS numbers: 71.10.Pm, 72.80.Sk 



In conventional BCS superconductors the amplitude 
of the complex order parameter |A|e l * vanishes at the 
transition temperature T c . This is in contrast to the un- 
derdoped cuprates, where experimental evidence (T), Q 
suggests that the transition is instead driven by the dis- 
ordering of the superconducting phase through thermal 
fluctuations, while leaving the magnitude |A| of the or- 
der parameter intact. A quantitative measure for the 
strength of phase fluctuations is provided by the ratio 
Q = 2T c /np s (0), where p s (0) is the zero temperature 
phase stiffness. This ratio determines how close the tran- 
sition is to being mean-ficld-likc. In BCS superconduc- 
tors Q -C 1, while in the underdoped cuprates Q ~ 1 
||. The effects of thermal phase fluctuations on d-wave 
superconductors have been investigated before, see e.g. 
Refs 13-flol. A key objective of these works is to identify 
clear signatures of thermal phase fluctuations in single 
particle properties such as the spectral function measured 
by ARPES and STS (Scanning Tunneling Spectroscopy). 
The purpose of this note is to provide an exact result for 
the partially integrated spectral function of a phase fluc- 
tuating superconductor in a particular limit. The latter 
quantity is defined as 



J dfc||A( W ,k) 



(1) 



where k\\ , k± are the wave vector components parallel 
and perpendicular to the Fermi velocity at the point of 
observation. We start by summarizing the essential as- 
sumptions underlying the model proposed by one of the 
authors and M. Khodas in 0, HH. The starting point 
is a superconductor with a general order parameter that 
arises from pairing on a Fermi surface the shape of which 
we keep general for now. In particular it could be open or 
consist of several pockets, as is believed to be the case in 
underdoped cuprates [F3 - tI5l |. Our following analysis is 
based the existence of well defined quasiparticles, which 
is a reasonable assumption for the nodal regions. The 
corresponding Bogoliubov-deGennes Hamiltonian is 

H = Jdvdv' [* + (r)] T {<5(r-r')e(-iV)r 3 

+ iA(r,r')r++h.c.}*(r'), (2) 



where we have defined Nambu spinors 1 i< T = (ipf, r a 
are Pauli matrices and the pairing amplitude can be cast 
in the form 



A(r,r') = A(r-r') e : < 



(R) 



(3) 



Here A(r) determines the symmetry of the order param- 
eter and R = (r + r')/2 is the centre of mass co-ordinate. 
Following the standard assumptions we neglect quantum 
fluctuations of A and focus exclusively on thermal fluc- 
tuations of the phase <p. The key point is to choose an 
appropriate model for these phase fluctuations. The ef- 
fects of fully three dimensional fluctuations are well stud- 
ied in the literature [r| and are found to be small. On 
the other hand, one would expect the spatial anisotropy 
of layered materials like the cuprates to strongly enhance 
the role of phase fluctuations. The extreme limit would 
be the purely two dimensional case, on which we fo- 
cus in what follows. We emphasize that even purely 
2D models have a window of applicability to e.g. thin 
films [l?! and x = 1/8 LBCO, where the phase transi- 
tion was found to be of Berczinskii-Kostcrlitz-Thouless 
(BKT) type , G! , Hil . Similarly, the analyses of the 
temperature dependence of magnetization, Lond on p en- 
etration [2l[ depth and terahertz conductivity j22l | for 
high quality underdoped BiSCO crystals show that al- 
though the superconductivity below T c is of a 3D nature, 
the superconducting transition in these systems is rather 
close to a BKT transition. In the latter case our the- 
ory will be applicable in a temperature regime above T c , 
where the phase correlation length is exponentially large 
and the phase fluctuations can effectively be considered 
as critical. 

In the mean field approximation fluctuations of the 
order parameter A arc ignored, and the resulting Green 
function takes the familiar BSC form 



G B cs(w,k) 



(w + iO) 2 -e 2 - A 2 (k) 



(4) 



The corresponding spectral function — -Im Gbcs con- 
sists of two delta function peaks centered at positive and 
negative frequencies. These peaks will be broadened by 
thermal phase fluctuations. The following facts are of 
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crucial importance in the following: (i) Since the long 
wavelength fluctuations are classical, the electron fre- 
quency is conserved, (ii) Since in the region of interest 
the amplitude |A| is assumed to be fixed, self-consistency 
between the electron Green function and the order pa- 
rameter is not an issue. Hence the calculation of the 
spectral function is reduced to solving the Bogoliubov- 
deGcnncs equations for a particle with pairing amplitude 
([3]) and then averaging the result over a given distribution 
of phase fluctuations, (iii) Since we are interested only 
in long wavelength fluctuations, the distribution function 
P{4>) = c~ F *l T can be fixed by symmetry considerations: 
as long as the discrete lattice symmetries include C4, the 
group of in-plane rotations by ninety degrees, the distri- 
bution function must be spatially isotropic (apart from 
irrelevant higher gradient terms). This leads to 



T 



2T 



dxdy [{d x <j)f + (d y ^f 



(5) 



where the prefactor T _1 results from the integration over 
imaginary time. In contrast to the phase fluctuations the 
Green's function at low energies is very different in the 
directions perpendicular and tangential to the Fermi sur- 
face. In a process where an electron close to the Fermi 
surface changes its momentum from k to k + q by scat- 
tering off the pairing potential its Green's function is 

Go^k + q) = w-e(k + q) 

« w -e(k)-«g||- |I. (6) 

Here Qy and q± are the components of the momentum 
respectively parallel and perpendicular to the Fermi ve- 
locity Ve(k). As a result of the isotropy of the distribu- 
tion of phase fluctuations the typical values of q\\ and 
q_i arc the same and of order A max (the maximal gap). 
Therefore the last term in ([5]) is proportional to the small 
parameter A/ep. If we neglect such small corrections the 
electron propagates along a straight line in real space 
and the transverse momentum is conserved. The elec- 
tron Green's function can then be calculated separately 
for each frequency w and Fermi surface point k. Under 
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FIG. 1: Geometry denning the model in l[7|. 
the assumptions summarized above the initial problem 



([2]) is recast as a field theory described by the Lagrangian 
L = Fh, + f,^5£ with 



-iu> n I - ivr z d x 



A(k x ,x) , , A*(k x ,x) 



where A(k±,x) = A(fc^)c 1< ^^' ' and we have intro- 
duced <3> Wri = (Vw.Ti v—un i) T ■ I n we have used 
co-ordinates as shown in FigfTJ As was pointed out in 
Ref.Q, the model defined through Eqs and ([5]) is in 
fact equivalent to the anisotropic spin- 1/2 Kondo prob- 
lem. In terms of this impurity model the phase fluctua- 
tions play the role of the host, while a single Bogoliubov 
quasiparticle constitutes the magnetic impurity. The re- 
duction of the underlying interacting electron model to 
a single-impurity problem is possible because the emer- 
gent low-energy degrees of freedom are non-interacting 
Bogoliubov quasiparticlcs. The many-body aspects of 
the problem enter the determination of |A|, but as this 
is treated as a parameter of our model we can avoid the 
issue of its calculation. Under a field redefinition 



X. 



(7) 

and subsequent analytic continuation \uj n — > uj + iO we 
obtain the Hamiltonian 



H cS = v h{Lu + i0)f 3 + H hum [4>] 



(8) 



A(fcx) 
2v 



where f a = x +7 " a X i s a short hand notation for fermionic 
bilinears. In this setting the coordinate x plays the role 
of Matsubara time. It is dual to the momentum com- 
ponent fc|| parallel to the Fermi velocity at the point of 
observation. We note that in the approximation under- 
lying ((7]) the electron momentum parallel to the Fermi 
surface is conserved so that the fcrmions \ depend only 
on x, while the phase field is a function of both x and 
y. For convenience we assign \ the coordinate y = 0. 
Since the fermion number is conserved, the f-operators 
are in fact components of a spin S=l/2. 

The Hamiltonian i?buik arising from ([5]) describes the 
phase fluctuations. For temperatures below the BKT 
transition temperature Tbkt == only smooth field 

configurations contribute so that 



1 f°° 

ffbulkfe] = J dy [(4nd) 2 U 2 



{d y <t>?\ 



(9) 



where n is the momentum density conjugate to the 
field 4>, with equal time commutator [n(yi), 4>(y2)\ _ = 
—\6{yi — 2/2)- In order to be able to treat the tempera- 
ture region T > Tbkt we need to allow singular (vortex) 
configurations of the (j> field. The effects of vortices can 
be illustrated for the example of the two point correla- 
tion function of bosonic exponents . The latter takes the 
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form 



( e l<«>-i) e -l<«r2)\ 



2d 



F L(T)) ' 



(10) 



where d = T/(8Tbkt) is the scaling dimension of the 
order parameter, £(T) is the correlation length and b ~ 
(v/ep) is the short distance cut-off. The short and long 
distance behaviour of the scaling function is F(p <1) = 
p~ 2d and F(p > 1) ~ K {p) respectively (see also [23| ) . 
Below the transition (where £ = oo) the function (jTUJ) 
decays as a power law and above the transition where the 
vortices are relevant it decays exponentially with finite 
correlation length £(T). We show below how to take this 
into account. It was shown in Ref [lij that jHJ) , ([9]) is 
equivalent to the anisotropic S=l/2 Kondo model in the 
regime of extreme anisotropy g» ^> g± 



(11) 



J 



■i + z , g± 



N 



a t* T aa> a P°' + h - C - 



p.k 



where the magnetic field h is related to the real frequency 
lu in ([8]) by analytic continuation h = iu> + 0. Our main 
result derives from the observation that the partial den- 
sity of states (PDOS) defined by |T]) is equal to the Green 
function of x fermions at coinciding coordinates x. Tak- 
ing into account the change of variables ([7]) we find that 
the PDOS is obtained by analytic continuation of the 
impurity magnetization of the Kondo model (jlip 



p P (uj)/po = 2Re M(h = iuj + 0) 



(12) 



where ira = 2 (i-d) M lnrf+(l — d) ln(l — d)]. As a function 
of a complex variable M(z) admits a power series expan- 
sion in odd powers of z for \z\ < 1 and concomitantly is 
purely imaginary along the imaginary axis. By virtue of 
the identification (jT2"j) this implies that the PDOS van- 
ishes at |w| < Th- Thus there is a sharp gap equal to Th 
in the density of states, which at T ^ is always smaller 
than the mean field gap A(k±). On the other hand, for 
\z\ > 1 the following expansion holds 



Po 



1 



sin[27rn<i] 
^ 27r 3 / 2 (n!) 

n=l v ' 

T H e- ,ra i 2n ( 1 - d ) 



E 



r(nd)r(i + (l-d)n) 



We note that 



lim 



Pp(uj) 



d^O po 



A 2 



A) , 



(15) 



(16) 



which corresponds to the mean field result. In order to 
establish the relation between the gap Tjj and A(k±), d 
we compare (|15[) to the perturbativc expansion for the 
PDOS in the model (JSj) . Second order perturbation the- 
ory gives 



Spp 
Po 



cos(7rd)r(2 - 2d)2 2d A 2 b 2d 



2w 2 (!- d ) 

which yields the desired identification 
T H = A(k ± )VT~d\VdbA(kj_)] d/1 d \r(l-d) 

(18) 

Given the result (fT5|) for the PDOS we may calculate the 



(17) 



l/l-d 



Here po is the bare density of states. This expression 
provides a link between spectral properties of the sin- 
gle electron problem ^ and thermodynamic properties 
of the many-body theory (jTTJ) . In order to utilize the 
known exact expression for M(h) in the Kondo problem 
[25| we need to relate the parameters Jgu, Jg± in (fTl) 
to d and A(fcj_). The interactions in the Kondo model 
increase under renormalization and enter the strong cou- 
pling regime at a scale Th which is known from the exact 
solution [23 



T 



eF(5±/.9||) 2V9 " 



(13) 



On the other hand the usual scaling argument gives Th ~ 
epg 1 ^ 1 d \ which leads to the identification g\\/2ir = 1—d 
with d = T/(8Tbkt)- The expression for the impurity 
magnetization derived in [25[ then reads (the parameter 
p in [25[ is related to d by /x = w(l — d)): 



M(h/T H ) = 



dx r i-i^a r l + Lr 



rfi-i^ 



4tt 3 / 2 J_ oq x + iO 
x exp|-2ix ]n(h/T H ) + na |, (14) 
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FIG. 2: Partial DOS as a function of frequency plotted for 
different temperatures T = 8d Tbkt- The ratio of A(fej_) to 
the cutoff 1/6 is fixed as 0.1. Due to particle-hole symmetry 
Pp{~uj) = pp(uj). 



full tunneling density of states, 
superconductor this results in 

p(uj) cx 



In the case of a d-wave 



(19) 
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In Fig. [2] we show the PDOS ([T5]) as a function of fre- 
quency for several different temperatures. The most no- 
ticeable feature is the persistence of a sharp gap. In ad- 
dition we observe that the singularity characteristic of 
the BCS mean- field solution is strongly suppressed as T 
increases. This demonstrates that thermal phase fluctu- 
ations have a sizable effect on integrated spectral proper- 
ties. In realistic materials the sharp gap will be smeared 
by both impurity scattering and the effects of Fermi sur- 
face curvature neglected in our analysis. In Fig. [3] we 




bA(k.) 



FIG. 3: The dimensionless gap Th b as a function of q = 
A(k ± )b for d = 0.05,0.125,0.25. The lowest curve, corre- 
sponding to T = 2Tbkt, has been modified to indicate the 
effects of vortices as described in the main text. 

plot the temperature evolution of the gap Tjj- We see 
that temperature effects are negligible. In particular this 
implies that the (d-wave) form of the gap remains robust 



through the transition, at least if vortices are ignored. As 
discussed above, the main effect of vortices is to induce 
a finite correlation length £(T). In the corresponding 
Kondo picture this translates to a gap in the excitation 
spectrum of the host. Exact results are available in this 
case as well [13] ■ On a qualitative level what happens in 
the Kondo picture is the following: as long as the corre- 
lation length £(T) is larger than the inverse Kondo scale 
v /Th the vortices have little effect on the physical prop- 
erties. However, as soon as £(T) falls below v/Tg the 
scaling terminates before the strong coupling regime is 
reached. As a consequence the gap in the PDOS is re- 
duced for momenta close to the node k± < 1/[6A£(T)]. 
We have indicated this effect in the dotted curve in Fig 13] 

In this work we have considered a model for thermal 
phase fluctuations in a superconductor recently proposed 
in Ref. 0] ■ By exploiting a mapping to an effective spin- 
1/2 Kondo problem we have derived an exact result for 
the partially integrated spectral function (JTJ). Our main 
result is that thermal fluctuations have a substantial ef- 
fect on the single particle spectral function. The best can- 
didate for comparing our theory to experiment is x=l/8 
do ped LBCO. ARPES and STS measurements performed 
in [26( show that the d-wave gap is already well formed 
at the BKT transition. It would be interesting to map 
out the detailed temperature dependence of the spectral 
function by ARPES in the region of strong diamagnctic 
fluctuations T < 40K and carry out a partial integration 
along the nodal direction. 
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